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Most stars will experience episodes of substantial mass loss at some point in 
their lives. For very massive stars, mass loss dominates their evolution, al- 
though the mass loss rates are not known exactly, particularly once the star 
has left the main sequence. Direct observations of the stellar winds of massive 
stars can give information on the current mass-loss rates, while studies of the 
ring nebulae and HI shells that surround many Wolf-Rayet (WR) and lumi- 
nous blue variable (LBV) stars provide information on the previous mass-loss 
history. The evolution of the most massive stars, M > 25M Q , essentially fol- 
lows the sequence O star — > LBV or red supergiant (RSG) — > WR star — > 
supernova. For stars of mass less than ~ 25Mq there is no final WR stage. 
During the main sequence and WR stages, the mass loss takes the form of 
highly supersonic stellar winds, which blow bubbles in the interstellar and cir- 
cumstellar medium. In this way, the mechanical luminosity of the stellar wind 
is converted into kinetic energy of the swept-up ambient material, which is 
important for the dynamics of the interstellar medium. In this article, analytic 
and numerical models are used to describe the hydrodynamics and energetics 
of wind-blown bubbles. A brief review of observations of bubbles is given, and 
the degree to which theory is supported by observations is discussed. 



1 Classical wind-blown bubbles 

The effect that a stellar wind has on its surroundings can be dramatic. A 
massive star will first have formed an H II region around itself, and so the 
stellar wind interacts with this ionized gas. A fast (typically 2000 km s _1 ) 
stellar wind is hypersonic with respect to the ambient medium (sound speed 
~10 km s _1 in photoionized gas) and so a two-shock flow pattern forms: one 
shock sweeps up the ambient medium, accelerating, compressing and heating 
it, while th e other shock decelerates the stellar wind i tself, heating and com- 
pressing it llPikel'neil Il968l iDvson fc de Vriesl Il972t lAvedisoval Il972ft . The 
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Fig. 1. Schematic of the different regions of a wind-blown bubble around a massive 
star. i? s i is the stellar wind shock, R s is the shock set up in the ambient medium, 
and _Rs3 is the isothermal shock sent out ahead of the ionization front that borders 
the H II region, while R c is the contact discontinuity separating shocked stellar 
wind material from swept-up ambient medium and Ri marks the ionization front, 
separating neutral and ionized gas. 



basic scenario is shown in Figure ^ I n most treatments, the photionized re- 
gion and the neutral shell beyond the stellar wind bubble are not discussed 
separately, and it is assumed that the shock R s expands into a fully ionized 
ambient medium. 

Conditions behind the ambient medium shock, i? s , favour rapid cooling 
of the swept-up shell of material, which occurs once the shock velocity fall s 
below 200 km s _1 , so that this region becomes thin and dense llFalleLll97.^ . 
The shocked stellar wind, however, can have temperature > 10 7 K and thus 
does not cool efficiently. A hot, very low density bubble of shocked stellar 
wind forms, which pushes the cold, swept-up shell like a piston. Since the hot 
shocked wind remains adiabatic, this is known as an "energy-driven" flow. 

If the ambient medium has a radial power-law density distribution p(r) — 
Por 13 , and R s and V s = R s , are the radius and velocity of the bubble-driven 
shock, respectively, then the expansion law for the hot bubble can be found 
by considering momentum and energy conservation, with the assumpti on that 
R s = R r , i .e., that the region of swept-up ambient medium is thin I Dvsorl 
Il984l [l989^) . The momentum and energy equations are 
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d /4vr , P \ , 

dt{T R °—i)= E "-^ R * PV ° (2) 

where P is the thermal pressure in the hot bubble, giving as solution 

(if \ 1/{5+0) 

R s = 1 /( 5 +« ( ^ j t 3 '^ , (3) 



with 



(3 + /3)(5 + /3) 3 



127r(ll + /3)(7 + 2/3) ' 

Here, the stellar wind mechanical luminosity, ^ = ^M W V W , is taken to be 
constant, and M w , V w are the stellar wind mass-loss rate and velocity, respec- 
tively. A time-varying stellar wi nd mechanical luminosity was considered by 



iGarcia-Segura fc Mac The case f3 = (i.e., constant density am- 

bient medium) returns the well-know n formulae of the thi n shell approxima- 
tion for stellar wind bubble evolution ijWeaver et alllT977j) . The case f3 = — 2 
would represent a wind bubble expanding into the density distribution left by 
the wind of a previous evolutionary stage of the star. 

"Momentum-driven" flows occur when the shocked wind cools in a time 
much less than the dynamical timescale. The swept-up shell of ambient ma- 
terial is driven by the momentum of the wind rather than by the pressure of 
the shocked wind. The expansion law can be found by considering momentum 
conservation. 

j t {j-Po R3 S +0V s)=M w V w . (6) 

For constant wind momen t um, M m V,„ , the radius and velocity of the outer 
shock are given by IjDvsonl 1 1984 Il989j) 

/ • \ l/(4+/3) 

R s = x 1/(4+/3) * 2/(4+/5) , (7) 




with 

x =(i±m±£>. (9 ) 

A flow can be defined as energy or momentum driven by the thermal be- 
haviour of the shocked stellar wind at a time, to, when the wind has swept 
up its own mass of ambient material. If the shocked wind can cool by this 
time then the flow will be momentum driven, otherwise it is energy driven. 
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The cooling time in the shocked stellar wind can be estimated using Kahn's 
approximation to the cooling rate, L = AT~ 1 / 2 n 2 erg cm~ 3 s -1 , in the tem- 
perature range 10 5 < T < 2 x 10 7 K (jKahnl Il976^ . where A = 1.3 x 1CT 19 is a 
constant and n and T are the number density and temperature, respectively, 
in the shocked gas. The cooling time is thus 

tcool = ~, ~~7T7" « 4 X 10- 35 -^ , (10) 
(7 - 1)£ PR 

where a mean nucleon mass of m p = 2 x 10~ 24 g has been used and it is 
assumed that the gas is fully ionized. 1 For a momentum-driven flow, pn = 
M w V w /4:irRs 1 , where i? s i is assumed to be the same as R s (see Eq. 0), since 
both the shocked wind and shocked swept-up ambient medium shells are thin 
in this approximation. Hence 

i coo] ~ l X lO-33 x 2/(4+/3) i ^-(2+/ 3 )/(4 + ^) F a8+^)/(4 + 0) p -2/(4 + /3) t 2 / (4 + ^ (n) 

in this case. 

When the flow is energy driven, the interior pressure of the hot bubble, Pb, 
is uniform and can be found by substituting equations[3]and0]into equation^ 
The density in the hot, shocked stellar wind is n = Pi,/2kT, where T is the 
post-shock temperature, which depends only on the stellar wind velocity. In 
this case, therefore, the cooling time is given by 

t cool ~ 2 X 10-35^-(2+«/(5 + /3)^3(7 + /3)/(5+/3) / ,-3/(5+/3) i (4-^)/( 5+/ J) ^ {u) 

where 

(5 + /3) 2 (3 + /3) 
? (7 + 2/3)0 2 +/3 ' 1 ' 

The ratio of the cooling time to the dynamical time thus varies as £~' 3 /( 4 +' 3 ) 
for a momentum-driven flow and t _ ( 1 + 2 ' 5 )/( 5 +/ 3 ) for an energy-driven flow. 
While the cooling time is greater than the dynamical time, an initially energy- 
driven flow can remain energy-driven. Thus, flows can remain energy driven if 
/3<— \. If > /3 > — h, energy-driven flows can become momentum driven. 
An initially momentum-driven flow can only be maintained if /3 > 0, otherwise 
it becomes energy driven. 

Main sequence stellar winds interacting with a uniform density ambient 
medium (/3 = 0) will normally be initially energy driven, unless the wind itself 
is very dense and can cool immediately once it shocks. The flow will remain 
energy driven because of the high postshock temperatures until, eventually, 
this gas begins to cool. If thermal conduction is important between the hot, 
shocked wind region and the cold swept-up shell, this can lower the tempera- 
ture and raise the density in the hot bubble. This would enhance the cooling 
rate and so the changeover to a momentum-driven flow would occur earlier. 



1 The strong shock conditions p = jPrV^, n = Apu/m v and T = ^V^rn p /k have 
been used in the derivation. 
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A Wolf-Rayet star wind has a high mass-loss rate and could initially pro- 
duce a momentum- driven flow. However, the density gradient of the red super- 
giant wind (j3 = —2) would quickly lead to a changeover to an energy-driven 
flow. 

Observationally, it should be possible to discriminate between energy- 
driven and momentum-driven flows by determining the efficiency factors for 
the conversion of stellar wind mechanic al energy and moment um into swept- 
up gas kinetic energy and momentum l|Treffers fc ChuL Il982l) . From the an- 
alytical model, an energy-driven flow has a kinetic energy efficiency factor 
e = 3(5 + /3)/(ll + /3)(7 + 2/3) and a momentum efficiency factor \x = e(V w /V s ), 
while for a momentum-driven flow the corresponding values are e = V s /V w 
and fj, — I. Values of e <S 0.1 are taken to indicate a momentum-driven flow. In 
practice, however, it is difficult to estimate e since th e masses of both neutra l 
and ionized components must be taken into account l|Dvson &: Smithl . ll985(l . 
Furthermore, outer shells may no longer be driven directly by the current 
stellar wind if the bubble has become depressurized. 



Thermal conduction 

IWeaver et al.l l)l977t) obtain a self-similar solution for the temperature struc- 
ture across the hot bubble under the assumptions that the pressure is constant 
with radius and thermal conduction from the hot shocked wind to the cold 
swept-up shell is important. It is also assumed that radiative losses are negli- 
gible compared to the conductive and mechanical energy fluxes. In this model, 
the mass within the hot bubble is dominated by the evaporated mass from 
the shell, while this mass is negligible compared with that remaining in the 
shell. The self-similar solution for the temperature distribution across the hot 
bubble is 

T(r)=T b (l--^) 2 / 5 , (14) 

where Tf, = aEw 35 n^ 35 t~ 6 / 35 is the temperature in the inner part of the 
bubble and a is a constant which depends on the value of the conduction co- 
efficient. This model has been extended to include the effect of expansion in 
a power-law density distribution and non-spherical (but uniform) expansion 
bv lGarcia-Seeura Ez Mac Low! l|l995aj) . but the same depe ndence on the sim- 
ilarity variable r/R s is retained. It has been argued fe.g.. iDvsori Il98lj) that 
thermal conduction will not, in fact, be important in stellar wind bubbles 
since the magnetic fields in the swept-up material will suppress it. It turns 
out that the temperature an d density profiles predicted by the self-similar 
model of IWeaver et alJ l)l977|) significantly overestimate the X-ray emission, 
as compare d to what has been o bserved, for a given set of stellar wind bubble 
parameters llWrigge et all 120051) . 



6 S. Jane Arthur 



2 Numerical simulations of wind-blown bubbles 

Numerical simulations are a powerful tool for studying the hydrodynamics 
of wind-blown bubbles, complementing analytical studies since many more 
physical processes can be included. 2 The first, spherically symmetric, numer- 
ical studies included radiative cooling and were able to follow the transition 
between the initial fully adiabatic s tage to the stage where the swept-up in- 
terstellar medium cools radiatively Subsequent numerical stud- 
ies in two dimensions show that instabilities develop during the evolution 
of wind-blown bubbles, where the nature of the instability present depends 
on the evolutionary stage and physics being considered. Kelvin-Hclmholtz, 
Ray leigh- Taylor, Richt meyer-Meshkov, Vishniac, and thermal instabilities 
have all been reported jR ozvczkal|l985tlRozvczka fc Tenorio-Taglel Il985albt 
Garcia-Segura fc Mac Lowul995lilGarcia-Segura et aJlll996ala : lBrighenti fc 
1995alblll997tlStrickland fc Stevenslli998HFrever et alJ.l200^r However. care 
must be taken in the identification of these instabilities, since the resolution 
of the numerical grid can play an important role, as can parameters such as 
the cutoff temperature below whi ch radiative cooling is assumed equal to zero 
( see IZhekov fc Mvasnikovl l2000bl for a more detailed discussion). Although 
the most widely used analytical models take into account thermal conduc- 
tion ([Weaver et allll97lllGarria-Segura fc Mac LowLll995al). very few numer- 
ical s t udies inc lude this process (notable exceptions are lZhekov fc Mvasnikovl 
Il998l l2000albh . The presence of even a weak magnetic field will be enough 
to cause asymmetric thermal conduction and will thus lead to an asymmetric 
wind-blown bubble l|Zhekov fc Mvasnikovl l2000b(l . 

In this section, the evolution of a wind-blown bubble around a 4OM star 
from the main-sequence phase, through the red supergiant phase to the final 
Wolf-Rayet phase is examined by means of numerical simulations. The first 
two stages are treated as one-dimensional (spherically symmetric) simulations, 
while the start of the Wolf-Rayet phase is modelled in two dimensions (cylin- 
dr ical symmetry ) . Pho toionization is included using the method described 



in iHennev et al. (2005). Radiative cooling and photoionization heating are 
also taken into account. Thermal conduction and magnetic fields are not in- 
cluded in these models. A simple three-wind model is considered, where the 
stellar wind mass-loss rate and velocity is taken to be constant during each 
phase and the ambient medium has uniform dens i ty and temperature (cf. 
iGarcia-Segura fc Mac Lowl Il995bt Ivan Marie et all 12005(1 . Although a par- 
ticular simulation is described here, the general features are common to all 
simulations of the evolution of wind-blown bubbles. 
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Fig. 2. Radial structure of a stellar wind bubble and H II region around a 40Mq star 
at an early stage of evolution (after 40,000 years). Top panel: total number density 
(solid line) and ionized number density (dotted line). Bottom panel: temperature 
(solid line) and internal energy (e = p/[y — 1] — dotted line). Stellar parameters 



9.1 x 10" 7 M, 



from Ivan Marie et al l 1200 



V w 



890 km s~\ and S, = 6.34 x 10 47 s _1 (adapted 



2.1 Main sequence evolution 

The main-sequence evolution of the ionized region (H II region) and stel- 
lar wind bubble around a 40M Q star in a uniform, neutral ambient medium 
with constant stellar wind and ionizing photon rate is presented here. The 
initial neutral ambient medium is taken to have density ?ih = 15 cm~ 3 
and temperature 300 K, while the stellar wind has a mass-loss rate of 
M w = 9.1 x 10~ 7 A/ Q yr _1 and wind terminal velocit y V,„ = 890 km s" 1 , 
and the ionizing photo n rate is S* — 6.34 x 10 47 s~ : (cf. ISchaller et all Il992t 
Ivan Marie et allhoOSh . 3 

The H II region forms immediately, reaching an initial Stromgren radius of 
~ 4.5 pc. A low-velocity, isothermal shock forms ahead of the ionization front 
in the neutral gas and the hot (T = 10 4 K) ionized region begins to expand. 



2 Numerical simulations solve the full gas-dynamic equations and physical processes 
are generally included via source terms. 

3 In these simulations, the stellar wind is injected in the form of a thermal energy 
and mass source in a small volume around the star, and for this reason tempera- 
tures in the unshocked stellar wind region are high but do not affect the dynamics 
since the flow is highly supersonic in this region. 
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The stellar wind bubble forms inside this ionized regi on. Cooling in the swept 
up ionized gas is very efficient (see, e.g.. iFallel Il975l) and this region quickly 
forms a thin, dense shell. The increase in opacity traps the ionizing photons 
and the main ionization front becomes confined to this swept-up shell. The 
region between the shell and the outer isothermal shock is occupied by gas 
that belonged to the initial Stromgren sphere and is now slowly recombin- 
ing (recombination time is approximately 10 5 /n years). This can be seen in 
Figure which shows the density, ionized density, temperature and thermal 
energy (i.e., pressure) structure after 40,000 years of evolution. The outer 
ionization front is no longer sharp since gas here has started to recombinc. 
However, the shock that was launched ahead of the ionization front continues 
to expand at < 10 km s _1 into the neutral ambient medium. 

In two-dimensional simulations, strong radiative cooling in the swept-up 
ambient medium leads to the formation of instabilities, wh i ch cau se corruga- 
tion of the shell l|Strickland fc Stevensl Il998l iFrever et all l2003|) . Shadowed 
regions then form in the H II region beyond the stellar wind shell and re- 
combination occurs faste r here, leading to a pattern of neutral sp okes and 
non-radial velocity helds l|Frever et all 1200:4 lArthur fc Hoard l2005h . 

Eventually, the wind bubble overtakes the outer neutral shock. By this 
time, the swept-up shell of material has spread out and its density has 
dropped. This is because the pressure in the shell is the same as that in the 
hot wind bubble, which falls slowly as the inner wind shock moves outwards. 
Since the temperature in the swept-up shell is determined by photoionization, 
and is thus a constant ~ 10 4 K, then the density has to fall (and hence the 
swept-up shell must broaden) to account for the fall in pressure. The inner 
wind shock will continue to move outwards until the ram pressure of the stel- 
lar wind here (p w V£ = M w V w /AttR 2 s ) balances the pressure of the interstellar 
medium. The H II region establishes itself in this, now, low-density shell, and 
sends a shock ahead into the neutral gas. 

Figure shows the radial structure of the stellar wind bubble and H II 
region at the end of the main sequence stage, after some 4.31 x 10 6 years of 
evolution. The inner stellar wind shock is located about 3 pc from the star. 
The hot (T ~ 10 7 K), very low density (m ~ 10~ 3 cm~ 3 ) shocked stellar 
wind bubble occupies the region 3 < R < 24.5 pc. Outside of this is the shell 
containing the H II region. The ionized gas has density rii ~ 1 cm -3 and 
temperature T = 10 4 K, and extends between 24.5 < R < 27.5 pc. The broad 
neutral shell, 27.5 < R < 33.5 pc, has density n n ~ 30 cm' 3 and temperature 
T ~ 10 3 K. The entire region between 3 < R < 33.5 pc has uniform pressure 
and is still expanding slowly into the ambient medium. 



2.2 Post-main-sequence evolution 

In Figurc lH the possible mass loss h istory of a 60M Q star is shown l|Langer et all 
ll994l:lG~arria-Segura et all Il996a|) . A detailed study of the evolution of the 
wind-blown bubble around stars of initial masses 60 M@ and 35 Mq has been 
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Fig. 3. Radial density and temperature structure of the stellar wind bubble around 
a 40Mq star at the end of the main sequence stage, i.e. after ~ 4.3 x 10 6 years. Line 
types are the same as for Fig. [5] 




2x10 s 
Time [yrs] 



3x1 l 



4x1 C 



Fig. 4. Variation of mass-loss rate (Mq yr and stellar wind velocity (km s x ) 
as a function of time for a 60Mq star (data from lGarcfa-Seeura et ailll996a|) . 
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carried out by Garcia-Segura et alJ lll996albl) . for the purely hydrodynamic 
case, and bv IFrever et alJ (|2003L l2005|) . including the radiative transfer of 
ionizing photons. The numerical simulations predict short-lived, observable 
nebulae during the LBV or RSG stage and the onset of the WR stage. Insta- 
bilities are formed in the dense shells swept-up by the different stellar wind 
stages, and these are consistent with clumps observed in ring nebulae around 
some Wolf-Rayet stars. The velocity of the LBV or RSG wind is found to 
play a key role in the detailed structure of the ring nebulae formed during 
this process. 



Red supergiant phase 
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Fig. 5. Radial structure of a stellar wind bubble around a 40 Mq star at the end of 
the red supergiant stage, i.e. after 4.5 x 10 B years. Top panel: total number density 
(solid line) and ionized number density (dotted line). Bottom panel: temperature 
(solid line) and internal energy (e = p/[y — 1] — dotted lino). 



During the red supergiant (RSG) phase, the ionizing photon luminosity of 
the 40 Mq star drops to 3.1 x 10 41 photons s _1 , and the stellar wind velocity 
drops to 15 km s" 1 whi le the mass- loss rate increases to 8.3 x 10 -5 M Q yr _1 
ijvan Marie et all [2005). This phase lasts only 2 x 10 5 years (an LBV phase 
would have an even higher mass-loss rate but a much shorter timescale, see 
Fig-EJ. The slow, dense RSG wind expands into the structure formed by the 
main sequence wind and a thin, dense shell of shocked RSG material forms 
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ahead of the freely expanding wind (see Figure [SJ. The high thermal pressure 
in the hot main sequence bubble causes backflow of low-density gas towards 
the star, which shocks against the RSG wind. The reduction in the ionizing 
photon rate leads to the disappearance of the H II region as the shell of ionized 
gas recombines. The hot bubble remains ionized since the recombination times 
are long in this hot, very low-density gas. The RSG wind itself is neutral. 

Wolf-Rayet phase 
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Fig. 6. Formation of Rayleigh- Taylor instabilities in the Wolf-Rayet wind shell just 
before the interaction with the red supergiant shell. The instabilities form as the 
WR wind accelerates into the r~ 2 density distribution left by the RSG wind. The 
temperature of the gas is shown 8000 years after the start of the Wolf-Rayet wind. 
Only one quadrant is shown, and the star is located at (0, 0). 

The Wolf-Rayet wind expands into the structure formed during the previ- 
ous two phases. During the Wolf-Rayet phase, the stellar wind becomes faster 
(V w = 2160 km s^ 1 ) and remains strong (M w = 4.1 x 10~ 5 M Q yrs -1 ) and 
the stellar ionizing photon rate also increases (S"* = 3.86 x 10 47 s _1 ). The 
transition between the RSG and WR winds will involve a period of wind ac- 
celeration into the r~ 2 density distribution left by the RSG wind, and this 
situation is Rayleigh- Taylor unstable. 4 The WR shell consists of swept-up 
RSG wind material and shocked WR wind. It becomes Rayleigh- Taylor un- 
stable at an early time in the simulation, well before the interaction with the 
RSG shell. Dense, cold clumps of material formed by the instability lag behind 

4 In the numerical simulations presented here, acceleration of the wind occurs ini- 
tially as a consequence of the way in which the stellar wind is included in the 
calculation. 
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the main WR shell but the shell itself is not disrupted, since the instabilities 
cease to grow once the acceleration stage is over (see Fig. EJ) ■ 
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Fig. 7. Interaction of the Wolf-Rayet wind with the dense shell from the red super- 
giant stage. Ray leigh- Taylor instabilities form as the WR wind accelerates into the 
r~ 2 density distribution left by the RSG wind. The instabilities become enhanced 
when the WR wind collides with the dense RSG shell and ruptures it. The panels 
show the temperature of the gas after times (a) 15,000 years, (b) 18,000 years (c) 
21,000 years and (d) 24,000 years after the start of the WR wind. 

Figure [7| shows the evolutionary sequence of the interaction of the WR 
shell with the RSG shell. The dense, cold clumps formed due to Rayleigh- 
Taylor instabilities in the WR shell act as obstacles around which the hot, 
shocked WR wind has to flow. Large amplitude ripples form in the RSG shell, 
through which the WR wind eventually breaks out into the low-density bubble 
formed by the main-sequence wind. The dense, cold clumps formed by the 
Rayleigh- Taylor instabilities move outwards more slowly. These clumps will 
be subject to hydrodynamic and photoablation but the numerical resolution 
of the two-dimensional simulation is not able to follow these processes. The 



Wind-Blown Bubbles 13 



mixing process produces gas with temperatures ~ 10 6 K, which should emit 
X-rays. 

During the red supergiant phase, 16.5 Mq of material are lost from the 
star. All of this material is swept into clumps and filaments by the Wolf- 
Rayet wind, which are potential mass-loading sources of the WR wind, to be 
discussed later. 




Fig. 8. Radial structure of a stellar wind bubble around a 40 Mq star at the end 
of the Wolf-Rayet phase, i.e. after a total time of 4.8 x 10 6 years. Top panel: total 
number density (solid line) and ionized number density (dotted line). Bottom panel: 
temperature (solid line) and internal energy (e = p/[j — 1] — dotted line). 



Figure |H1 shows the final radial structure of the wind bubble around the 
40Mq star just before it explodes as a supernova, from a one-dimensional 
(spherical symmetry) numerical simulation. The low-density bubble has been 
repressurized by the hot, shocked WR wind material, a new shock is being 
driven into the neutral shell, and the H II region is starting to reform. The 
one-dimensional simulation cannot show the fate of the clumps and filaments 
formed during the initial interaction between the WR wind and the RSG shell. 

2.3 Energy evolution in a wind-blown bubble 

Figure|H|shows the evolution of the neutral and ionized gas kinetic and thermal 
energy fractions with respect to the total stellar wind mechanical energy, E w t, 
resulting from the numerical simulations described above. During the main 
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Fig. 9. Top panel: Kinetic and thermal energies of the ionized and neutral gas as 
fractions of the total stellar wind mechanical energy as a function of time for the 
numerical simulation described in this paper. Thick lines represent kinetic energy 
and thin lines represent thermal energy, while solid lines represent neutral gas and 
dotted lines represent ionized gas energy fractions. The onset of the RSG and WR 
are indicated. Bottom panel: Kinetic and thermal energies as fractions of the WR 
stellar wind mechanical energy, E w (t — twrt) during the WR phase. 



sequence stage, the initial 400,000 yrs show large variations in the relative 
importance of the kinetic energies in the neutral and ionized gas. This is 
because this period of time sees the formation and expansion of the initial 
H II region, which then retreats back into the shell of material swept up by 
the stellar wind before expanding outwards again driven by the pressure in the 
hot bubble. The majority of the main-sequence stage sees an energy-driven 
bubble in which the thermal energy of the hot bubble pushes the expansion of 
a neutral shell of material. The kinetic energy of the swept-up neutral shell is 
roughly 0.25 of the stellar wind mechanical energy for most of this stage. This 
is la rger than the fract i on ^| predicted by the simple thin-shell model (see 
e.g.. lDvson fc Williamsl Il997j) and can be attributed to the H II region. It is 
the H II region that results in the different fractions summing to more than 
unity, since the contribution to the thermal energy by the heating of the H II 
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region in the shell surrounding the low-density bubble is substantial (see also 
iFrever et all l2003j) . Observations of only the ionized gas during this phase 
would show e <C 0.1, suggesting a momentum-driven flow, while observations 
of the neutral shell would give e > 0.1, hence an energy-driven flow. 

When the red supergiant phase begins, the ionizing photons are trapped 
close to the star and the outer H II region recombines. Also, the driving 
pressure in the hot bubble drops because the RSG stellar wind is such low 
velocity. As a result, both the thermal and kinetic energy fractions drop. 
Observations of both neutral and ionized gas during this phase would give 
e < 0.1, suggesting a momentum-driven flow. 

During the Wolf-Rayet phase, the RSG wind material is ionized by the high 
flux of ionizing photons and the high velocity WR wind forms a hot, high- 
pressure bubble. This is reflected in the energy fractions of the ionized gas, 
which rise steeply during this phase. It is a long time before the wind-blown 
bubble becomes pressurized out to the neutral shell, and so the neutral gas 
thermal and kinetic energy fractions remain low. Observations of the ionized 
gas during this stage would show e > 0.1 initially, but later in the evolution 
e <C 0.1 even though the flow is now energy driven. The neutral gas kinetic 
energy fraction corresponds to the outer swept-up shell, which is not repres- 
surized or accelerated during the timescale we are considering. Consequently, 
observations of the neutral shell during this stage of evolution would indicate 
that e < 0.1. 

FigureElshows that neither observations of only the ionized gas or only the 
neutral gas can reveal the full picture of the energetics of a wind-blown bubble. 
Observations of the outer neutral shell during this stage do not provide a true 
picture of what is happening with the interaction between the current stellar 
wind and its surroundings. 

2.4 Mass loading 

The presence of dense clumps embedded in a fast flowing stellar wind and 
subject to a large flux of ionizing photons from the central star suggests the 
possibility that material from the clumps will be incorporated into the flow 
by one or other mass-loading processes, namely photoevaporation, conductive 
evaporation or hydrodynamic ablation. If the amount of mass incorporated 
into the flow becomes important then the dynamics of the flow can be radically 
altered. Both supersonic and subsonic flows tend to an average transonic flow, 



if the mass loading is strong enough, the density in the hot shocked stellar 
wind can increase sufficiently that cooling becomes important in this region. 

Mass loading by embedded clumps in a stellar wind was first proposed to 
explain the ioniz ation potential-velocity correlation observe d in the WR ring 
nebula RCW 58 l|Smith et all ll984; Ha rtauist et allll986j) . In this scenario, 
the mass-loading is due to hydrodynamic ablation and has two regimes: one 




Furthermore, 
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in which the flow around the clump is subsonic and the other when it is su- 
personic. In the subsonic case, the mass-loading rate is proportional to M 4//3 , 
where M is the flow Mach number. In the supersonic case, the mass-loading 
rate saturates and is taken to be constant. Numerical modelling of the mass- 
loaded flow in RCW 58 shows that the position of the stellar wind shock is 
fixed by the mass lo ading and that coo l ing o ccurs in the outer parts of the 
hot shocked bubble ij Arthur et all Il99.4 Il99filb 



3 Observations of wind-blown bubbles 

Bubbles around massive s tars can be detected in a variety of ways: as optical 
ring-shaped nebula e fe.g.. IChul 1981). as she ll-shaped thermal radio contin- 
uum sources (e.g., iGoss fc Lozinskaval Il995|b as neutral gas voids and ex- 
pand ing shells in the hydrog en 21 cm line-em ission distribution ( Cappa et all 



2005), as infrared shells ( e . g. . IM a rstonl 1 1 9961) . and in a limited number of cases 
as diffuse X-ray sources (e.g., IWriggel 11999^1 . Gamma-ray burst afterglows 
probe the innermost regions of the progenitor wind bubble l|Chevalier et all 
|2004|) . and Type II and Type Ib/c supernova narrow spectral absorption fea- 
tures can be attribu ted to moving shells around the massive star progenitor 
l|Dopita et all Il984^ . 

Multiple concen tric shells of ma terial have been observed around many 
galactic WR stars l|Marstonl ^396). 5 The outermost shells are thought to 
correspond to the O star phase, and are most readily observable at far infrared 
or 21 cm radio wavelengths, although around 8% are seen optically. These 
outer shells have can very large diameters (> 100 pc), depending on the 
ambient ISM d ensity, and are expanding slowly (generally, v cxp < 10 km s _1 
IMarstonl Il996|) . The innermost shells are easily seen in optical narrowband 
images, and represent dense material recently ejected from the star in a LBV 
or RSG phase or swept up by the stellar wind in the WR phase. 

Of the ~ 150 observ ed Galactic WR stars about 25% are associated with 
optical ring-like nebu l ae faeckathorn et al.lll983lChu et allll983tlMiller fc Chul 
Il993l iMarston et all Il994alrjlb However, only ~ 10 of these WR ring nebu- 
lae have the sharp rims and short dynamical ages (ring radius divided by 
expansion velocity) that suggest that they are bubbles formed by the winds 
of the central stars during the WR phase. 6 The remaining ring nebulae are 



Ring n ebulae have also bee n observed around WR stars in the Large Magellanic 
Cloud iDopita et al.llT994l) but their distance makes them less easy to study. 
These objects are: S 308, RCW 58, RCW 104, NGC 2359, NGC 3199, NGC 6888, 
Anon (MR26), G2.4+1.4, and there is kinematical evidence that the nebulae 
around WR 11 6 and WR 133 have exp ansion velocities consistent with wind- 
blown bubbles (lEsteban fc Rosadol Il995ft. AnonfWR 128 ) and Anon(WR 134) 
are possibly wind-blown bubbles, too lGruendl et ai!l200ol) . There are also wind- 
blown bubbles around two other evolved stars: NGC 6164-5 (around an 06.5fp 
star), NGC 7635 (around an 06.5111 star). 
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less well defined and have dynamical ages much larger than the lifetime of a 
WR phase, and hence are thought to be simply stellar ejecta from the RSG 
or LBV stage photoionized by the WR star (like a planetary nebula), rather 
than dynamically sha ped by its w ind. Wind-blown bubb les were established 
as a class of object bv IChul l|l98lj) and lLozinskaval JT982). 
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Fig. 10. (a) Optical Ha and (b) [O III] A500 7 images of the ring n ebula RCW 58, 
together with (c) an image of their difference IIGruendl et all 12000). 



A striking example of an optical WR ring nebula is RCW 58 (see Fig- 
ure This nebula consists of clumps and filaments extending radially out- 
ward from the cent ral WN8 star, seen in Ha, w i th mo re diffuse [O III] emission 
extending beyond l)Chul Il982l iGruendl et all l20fl(ll) . The nebula is enriched 
in He and N, which indicates that the material has been ejected from the star, 
probably during a red supergiant stage. The expansion velocit y of the shell 
is 87 km s _1 , and broad linewidths indicate supersonic motions l|Smith et all 
Il988l) . The measured electron densities of the clumps are ~500 cm -3 . 

Most WR wind-blown bubbles are aspherical to a lesser (e.g., RCW 58, 
S 308, NGC 6888) or greater (e.g., NGC 2359, NGC 3199, G2.4+1.4) degree. 
The asphericity is often more pronounced in Ha images than in [O III], sug- 
gesting that it arises in the RSG or LBV ejecta phase. Such morphologies could 
also be due to the environment: for example, NGC 2359 and NGC 3199 are 
each bounded on one side by molecular clouds (|Sr.rmer)S et allh QSltlMaVstori 
l200ll) and their roughly parabolic shapes in Ha images could be the result 
of the wind-bubble blowing out from near the surface of the dense cloud and 
expanding d own a strong density gradien t, as has been suggested in the case 
of G2.4+1.4 IjDopita fc Lozinskavaill99(l) . An alternative explanation for this 
sort of morphology could be that the star is moving supersonically through 
the ambient medium and forms a bowshock ahead of it (this w as suggested 
as an explanation for G2.4+1.4 bv lBrighenti fc D'ErcoleU l995b). 

Four of the 8 WR wind-blown bubbles have been observed with X-ray 
satellites but only two of these, NGC 6888 and S 308, have been detected in 
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Fig. 11. Optical [O III] A5007 imag e (left) and XMM -Newton EPIC image (right) 
of the northwest quadrant of S 308 ijChu et alJ l)2003l ) : figure reproduced with per- 
mission from Revista Mexicana de Astrononomia y Astrofisica). 



diffuse X-rays. Both were detected with the ROSAT PSPC in the energy range 
0.1-2.4 keV and spectral analysis shows that the hot gas in these bubbles is 
dominated by the co mponent at ~ 1.5 xlO 6 K I Bochlrarevlll988tlWrigge et all 
ll994tlWriggelll999(l . although this could be partially contaminated by higher 
energy point source emission. ASCA SIS observations of NGC 6888 suggests 
that there is an addi t ional component of hot gas at a temperature of 8 x 10 6 K 
jWrigge et all Il99ct 1200^1 . while XMM-Newton EPIC observations of S 308 
indicate t hat the spectrum i s very soft, suggesting plasma temperatures of only 
1 x 10 6 K (|Chu et alll2003l) . Figures ITD and IP21 show the o ptical [O III ] A500 7 
and X-ray images of S 308 and NGC 6888, taken from IChu et~afl l|2003|) . 
The S 308 images clearly show that the X-ray emission is limb-brightened, 
and NGC 6888 appears to be limb-brightened also. 7 In both cases the X-ray 
emission is interior to the [O III] emission, which marks the position of the 
main shock wave. 



4 Observations confront theory 

It has been proposed that nebulae around massive stars begin with wind- 
blown shells in the main-sequence stage, evolve into amorphous H II regions, 
then ringlike H II regions, then nebulae with stel lar ejecta in the RS G or 
LBV phase before the star explodes as a supernova l|McKee et allll984|) . The 
bubbles blown by WR winds occur within the circumstellar material of the 

7 Although not clear from the ROSAT image, this is evident in recent Chandra 
images: http: //www. chandra . harvard . edu/photo/2003/ngc6888 . 



Wind-Blown Bubbles 19 



I | I I .1 



' • [out]. 



24 



1 



38 16 




y i 1 i 1 t i 1— 



' " ' 'J 

ROSAT PSPC 



4 



12 m 

Right Ascension (J2000) 



12 m 

Right Ascension (J2000) 



Fig. 12. Optical [O III1 A5 007 image (left) and ROSAT PSPC image (right) of 
NGC 6888 llChu et all 2003, figure reproduced with permission from Revista Mex- 
icana de Astronomia y Astrofisica) . 



RSG/LBV phase. In practise, wind-blown bubbles have rarely been observed 
around main-sequence O stars (although kinem atic eviden c e sug gests that 
wind-blown bubbles do exist inside H II regions iNaze et "all . l200ll) . Ringlike 
H II regions may be too faint to observe but expanding HI shells with masses 
up to 1O 4 M should occur around massive stars, as illustrated in Figure|21 In 
the later stages of massive star evolution and, indeed, around som e supernova 
remnants, expanding HI shells are detected l|Cappa et all l2005|) that must 
have originated in the main-sequence stage. 

The fact that so few of the known WR stars are surround ed by wind- blown 
bubbles suggests that this is a short-lived phase of evolution llChu lll98ll) . This 
is borne out by numerical simulations, which predict observable bubbles only 
during the early interacti on of the WR wind with the LBV or RSG nebula 
I Frever et, all Eooii . 120051) . More common are amorphous nebulae, which are 
observed around many WR stars, indicating a photoionization rather than 
wind-blown origin, even when there are clear indications that the central star 
possesses a strong stellar wind. These cases could be the result of the envi- 
ronment in which the nebula is found: if a strong density gradient is present 
the stellar wind will blow out in the direction of decreasing density ra ther 
than form a pressure-driven expanding bubble (|Arthur fc Hoard 20051 dis- 
cuss similar models in the context of cometary H II regions) . Clumps and fila- 
ments observed in nebulae such as RCW 58 are formed naturally in numerical 
simulati ons by Rayleigh- Taylor and other instabilities l|Garcfa-Segura et all 
1996a b). 

Theoretical estimates of X-ray emission as predicted by analytic theory 
fail to match the reality. Models that do not include thermal conduction 
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have very high temperatures in the shocked stellar wind and do not pro- 
duce the soft X-ray spectrum that has been det ected in NGC 6888 and 
S 308. On the other hand, the self- similar models of IWeaver et al 1 lll977h and 
iGarcia-Segura fc Mac Low! ljl995a|k which include thermal conduction, pre- 
dict much higher values of the soft X-ray luminosities than are observed. It 
appears that some sort of dissipative process , though not the classical thermal 
conduction proposed bv lWeaver et alJ l|l977h . is necessary to produce the soft 
X-rays. This argument presupposes that the X-ray emission is produced in the 
shocked Wolf-Rayet wind material. Two-dimensio nal numerical simula tions of 
the wind-blown bubble produced by a 35M Q star lFrever et alJ l|2005|) . on the 
other hand, show that the X-ray emission actually arises in the shocked, dense 
RSG wind material. This can be appreciated in Figured where gas tempera- 
tures of ~ 10 6 K, characteristic of soft X-rays are present in the shocked RSG 
shell material. Also, it could be expected that material photoevaporated or 
ablated from the dense clumps formed by the instabilities could augment the 
soft X-ray luminosity in the interaction region. 

5 Final Remarks 

Although the general characteristics of the formation and evolution of wind- 
blown bubbles around evolved stars are well understood, there remain a num- 
ber of important features to explain. Most obviously, it is interesting to know 
whether the asphericity of nebulae such as NGC 6888 and RCW 58 is due to 
the non-spherical ejection of material in the RSG or LBV stage, or whether 
it is due to some direction-dependent physical process, such as thermal con- 
duction, in the WR bubble phase. Photoionization is clearly important for 
the thermal balance and energetics of wind-blown bubbles and should not be 
neglected in numerical calculations. The role of instabilities in producing the 
clumps and filaments observed in some (but not all) bubbles should be exam- 
ined in more detail. Cooling within the clumps can lead to the formation of 
very dense, neutral structures within the hot, ionized bubble. The ablation, 
either by photoevaporation, thermal conduction or hydrodynamic ablation, of 
this neutral material will modify the physical properties and chemical abun- 
dances of the hot, shocked wind and such processes should, therefore, be stud- 
ied further. Theories should be tested with real data and new observations 
will always be needed, at all wavelengths, to provide more detailed kinematical 
information, abundance data, and improved stellar wind parameters. 
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